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CALCULATION OP ELECTROSTATIC FIELDS IN PERIODIC 
STRUCTURES OP COMPLEX SHAPE 

V. P. Kravchenko, V. I. Polevoy, A. M. Andrusenko 

and A. V. Kobyakova 


ABSTRACT. A method is presented for developing 
a universal algorithm for calculating electro- 
static fields in an infinite periodic structure 
consisting of electrodes of arbitrary shape which 
are located in mirror-symmetrical manner along 
the axis of electron-beam propagation. The 
method is based on the theory of R-functlons, and 
the differential operators which are derived on 
the basis of the functions. Numerical results 
are presented and the accuracy of the results is 
examined. The algorithm can be applied to the 
investigation of such devices as electron-optic 
systems, microwave-oscillator components, mass 
spectrometers, and linear accelerators. 

A method for developing a universal algorithm for 
calculating electrostatic fields in an infinite 
periodic structure consisting of electrodes of 
arbitrary shape which are randomly located in 
mirror-symmetric manner along the electron-beam 
propagation axis. The method is based on appll- 
catlon'-of the theory of R-functions. The accuracy 
of the results is evaluated. 


The problem of the distribution of electrostatic field in struc- 
tures with electrodes of complex form is of great practical value in 
electron optics. In the majority of cases, the complex geometry of 

* 
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electronic lenses does not permit the use of exact methods of 
solution; thus It becomes necessary to em lploy me thods of approximation. 
Since the numerical realization of these methods is usually accom- 
plished by means of an electronic computer, the choice of the method 
to be used should be made primarily in terms of its practical 
effectiveness. 

In the present paper, a method is proposed for developing a 
universal algorithm for calculating electrostatic fields in an 
infinite periodic structure consisting . of electrodes whose shape Is 
arbitrary and which are located at random In mirror-symmetric fashion 
along the propagation axis of the electron beam. 

The method is based on the use of special properties of R- 
functlons and the Introduction df the differential operators which 
are derived from them [1, 2]. R-functions are widely used in the 
solution of a broad range of problems in electrostatics and electro- 
dyanmics [3 - 7]- The merits of the present method, which is based 
on the above principle, lies in the possibility of describing analy- 
tically the shape of any plecewlse-smooth region and the broad 
applicability of the method to an entire class of similarly formu- 
lated problems. Moreover, the present method, as distinguished 
from the method of finite differences, requires a smaller expenditure 
of machine time and furnishes a solution in analytic form, which 
materially facilitates the analysis of the focusjing pr,bpertles of the 
structures under consideration. 

Given the following problem: to find the potential distribution 
of the electrostatic field in a certain structure consisting of an 
infinite set of electrodes, piecewise smooth in shape, located along 
the z-axls so that z„ — ni (u = 0, ±1, dc2...) are the planes of symmetry 
for the entire structure as well as for those electrodes which are 
intersected by these planes (Figure la, b, c). If inlthls case the 
potentials of the electrodes have the same symmetric form, then 
the condition: 

^ =0, (/i = 0,±l.±2...). (1) 

d\ ’n 


2 





Figure 1 


must be satisfied in the planes, where U is the desired potential 
of the electrostatic field, and v is the normal to the corresponding 
planes. 

The boundary condition (1) makes it possible to restrict the 
region to be considered along the z-axis to the planes and 
and to reduce the boundary problem to be considered to the equivalent 
problem for one "cell" of this structure of length I along the z-axls. 


Then the-'problem can be formulated as follows: in the region D, 

which is bounded with respect to z, the solution of the equation: 

&.U=0 (2) 

subject to the boundary conditions 

= 0, {i~\, 2, 3..., s), 

'%+» ^ ^ ^ 

where r< are the surfaces or boundaries of the electrodes for which 
the corresponding values of the potentials V. are given; T,-i-i is the 
surface or the boundary formed by the planes (axes) of symmetry; v 
is the internal normal to r,+» . 
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In those cases when the region D Is bounded (Figure la), or can 
be bounded (Figure lb) in the plane XOY, the solution for (2.) with 
the conditions (3) can be represented in the following structural 
form: 

b' = <l)„ + (o(D — (•!)„ + o.tl)), (h) 

where 


^0 ~ 


v k, n CO, 

f = l 

^ n o)/+ n w/ 

i+i 


( 5 ) 


<■»= 


•I 


fill ^^**>0 ^ I * fhtto 3 

• Om Ox ' Oy uy 


«?«llB 3 
Oz Oz ' 


where $ is a function of class subject to determination; w. is the /5 
function which gives an analytic description of the surface 1^^ of the 
corresponding electrode; Wg is the function which gives the analytic 
description of the surface formed by the planes of symmetry bound=:: 

Ing the region D along the z axis. 


It is not difficult to show that a solution in the form (^) 
subject to (5) satisfies the conditions (3), regardless of the form 
of the function;^.. For this purpose, it is sufficient to take account 
of the following.propertles of m., Ug, and the 'y5er^tor [1]: 


<o,|r. = 0; «)o|r,+,=0; % 

dv 


= 1 (/=1. 2.... s). 


(i>/ (M)>0; 0)0 (M)>0; 


( 6 ) 


DT (/)|r.+,= 


•1 

3\ 


*+l 


; D}«' (o)of) |r,+,=/jr,+,, 


■where M i s an arbitrary point of the region D, and f is any differ- 
entiable function. 


In the case the region D is not bounded in the plane XOY and the 
dimensions of the electrodes in this plane are finite, it becomes 
necessary to take into account the behavior of the field at infinity. 
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In this case, a more general form of the solution must be selected: 


U-<Do-r ■q '■ ~ u;'+ ir„ (o’ q ' 


(7) 


where 


W^o--T^: Q= }/ l+o; Qo= y l+wo: 


Qfn+l ■ " 

and m and mg are the degrees of the functions w and Wg, respectively, 


The solution In the form (7) not only satisfies the conditions 
(3), but also ensures the regularity of the function U at Infinity 
for any choice of the function fnom the class C^. 

The arbitrariness which Is permissible In the choice of the 
functi on"^ make s It possible to seek It by any method of approximation 
such that the structure (4) or (7) provides the best solution to 
equation ( 2) . 


Let us represent $ by means of the following approximation: 

p 

where ^ certain complete system of functions represented In 

the form: 


y.M = Pk (x) Pt(y)P<i (^) for the structure (4); 

= (f/ Pq (-) for the structure (7). 

On substutltlng (8) Into (4) and (7), and by virtue of the linearity 

of the operator we obtain: 

0 

U = Vo=- V • CwaY*,,. (9) 

*+i+*-o 

where, for the structure (4), 

(a,x*«). (10) 
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Analogously, we have for the structure (7): 

(OTto,). 


(10a) 


Next, we Introduce the function: 

p 

( 11 ) 

The problem (2) with the mixed conditions (3) can be reduced 
relative to V to the following problem with uniform boundary 
conditions : 

. = 0 . ( 12 ) 

'‘s+l 


AK=AVo: V'lr, = 0; 

dv. 


We shall seek the solution on the lineal of the functions V 

2 

belonging to C and satisfying the conditions (12). The given 
problem is equivalent to the problem of constructing an element of 
the energy space which realizes the minimum of the functional [8]: 

/ ( V) = (- AV. V) — 2{V, A'Fo) 5 (AKV+21'A'Fo) dx. ^ ^ 

On applying Green's theorem to (13) we have: 




dV aVo ■ (?V <3^0 , 
dx dx dy dy 




ClU), 


Taking Into account the conditions (6) and also for (7) the 
condition for the decrease of the po.tential at infinity, it is not 
difficult to verify that the surface integral In (l4) equals zero. 


Using the variational method of Rletz to minimize the functional ^ 
(1^), we obtain the following system of equations: 
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( 15 ) 


mmt ) \ Ox ux utj Oy 

*+»+<»“» « . 


Oi Oz / 


dx = 



dx ' Oy 


Oy 02 


dx. 

02 J 


Realization of the system (15) on an electronic computer permits 
the determination of the coefficient C*„ in the expansion (8), and 
therefore the solution of the problem posed. 


The proposed method is Illustrated below in an example of con- 
crete structures which are used in electron optics. 


Let us consider the structure (Figure la) consisting of two 
rows of bars with rectangular cross-section, infinitely long along 
the y-axis and distributed so that the planes Za=nl, {l=as) and the 
plane x = 0 are the symmetry planes of this structure. Moreover, 
above and below, the structure is bounded by the two conducting planes 
x=±bi . 

In particular, when the dimensions of the bars and the distance 
between them are the same 1^4=0 and V'i = K3 = — V2 , we obtain the struc- 
ture which is used for focusing the electron, current in generators 
of diffraction radiation [9]. 


Taking into consideration the homogeneity of the structure 
along the y-axls and its boundedness with respect to the x-axis, 
we shall seek the solution in the form (4) for a two-dimensional 
region D bounded by the lines 5-=0. z=a;, a- = 0, x = 6; (the solid lines 
in Figure la). 


For the given case: 




V' ;|i|;<ilV»4-<- V'7O)|>l)j0)«+ V'3ll)|<l>;0)<+ . 

+ <U|<1)2U>4 + UJ 1 W 2 WJ + W ’ 


( 16 ) 
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( 17 ) 


0) = (|)|<1)2(03104; «t»| = [/1A0M3]; 

j <1)2 = — c.)3 = — f/sAo/g/io]: 
i = — 2^0; jc = x — 63^0; fii = bj — 

/2 = x — 6)^0; /7 = ft4 — x^O; /, 2=2^0; 

/3 = &2 — x^O; /g=2 — a<>0; /i3 = a5 — 2>0; 

/4 = 2 — 02^0; /9=x — fes^O; fn = x^O, 

/5 = fl 3 — 2 > 0 ; /io = fr 6 — x> 0 ; 

Ag Is the symbol for the operation of R-conj unction [1] 

It can be verified by direct checking that (^) together with 
(16) and- (17) satisfies (3) for any choice of 

Representing $ in the form (8) arid using the method described 
above, we have for a system analogous to (15) where 'i'g and 
have the form (10). 


Let us consider a periodic structure (Figure lb) consisting of 
infinite plates of thickness 2d, the potentials of which alternately 
equal and l'^ plates there are openings formed by the 

intersection of these plates .with the surfaces of the form: 


_2r 21 

u-iT 






(r — ?«i )' 

.it 


V 


= 1 . 


(18) 


where zj, =n/ (n = 0. ±1, ±2...); / is the distance between the center of 
two neighboring plates, t is the exponent a. b. c, d are constants 

which determine the magnitude and shape of the openings, and c:S>d . 


/8 


The appearance of -an opening in (I8) is selected in connection 
with the fact that, depending on the magnitude of the exponent and 
the relationships among the quantities a, b, c, and d, the possibility 
emerges of obtaining openings of the most varied shape. Thus, if 
, we obtain an opening of practically rectangular shape (slit 
diaphragm). Then t = 1 and a = b, we obtain various types of 
axially dymmetric openings from segmental (when c = a = b) to cylin- 
drical (when c^^a. b. and d). 
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Consequently, from (l8) It is possible to obtain In particular 
cases the basic forms of diaphragms widely used In electrostatic 
focusing system. 


We shall seek a solution in the region D bounded by the planes 
z = 0, 2 =/, Jc=0, //=0 (the solid lines in Figure lb), given on them 
conditions of the type (1). Note that, beginning with some , 

the field between the plates may be considered as practically homo- 
geneous with a linear distribution of the potential along the z-axls. 
This makes it possible to restrict the region D In the plane XOY to 
a cylindrical surface of the form: 



(19) 


along z, where JCo^a, yo>b . 


Here on the surface (19) in the Interval d^z^(l — d) ', it Is ' . ' 
necessary to assign a distribution of the potential analogous to that 
which exists between two infinite planes. Taking this into account, 
the boundary conditions for the region D can be written In the 
following form: 


t/|r. = K.; iy|r. = l'2; 


I — 2d 



( 20 ) 


where F^, are the surfaces of the electrons; F^ is the surface 
(19) bounding the region D In the plane XOY, F|^ is the surface formed 
by the planes of symmetry. We shall select the structure of the 
solution In the form (4), where for the given case: 


<Do — 


V , 

0) j0l5 +* + 01 1 0)J + 10 


(I) — (|)|(|)2&)3: 0 )|= — [/1A0/2/3I; W2== [/lAoWs]; 

0)3 = 1/2/4 Ao/b] ; 0)0= [/t/sAo/oAo/ lol 

_ 2 / „2I , 2 < 

/|= ^ + jT “ f2 = d — z^0; 


a" r 

fs=d+z^0\ ft=z — l+d^{)-, ii==i+u — z^o. 


( 21 ) 
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x” I#" 

~5?H — 5 — 1^0; /7=2^0; 

*0 yl 

, /s= y (/ — 2)^0; /fl=x^0; /io=y^O; 

_ »'./.+ V'jf, 

^ -~2d-l ■ 

R epresenti ng <!>' in the form (8) and performing the sequence of 
operations described above, we obtain the system (15) to determine 
Ckig . The limits of integration with respect to z are from 0 to 1; 
with respect to x, from 0 to Xq; with respect to y, from 0 to yQ. 

It should be noted that the values Xq and y^ are selected in 
the process of numerical realization on an electronic computer, based 
on the condition that their subsequent increase must not influence 
(to a given degree of accuracy) the value of the,potentlal of the 
field U at any point ^oT the domain D under consideration. 

Let us consider the periodic structure (Figure Ic) consisting of 
a system of cylindrical rings with internal radii r^ and V 2 > thick- 
ness h^ and h 2 , length 2d^ and 2 d 2 , located along the z-axis at equal 
distances I between their centers. In the present instance, the 
region can be bounded only along the z-axls. Here the planes of sym- 
metry have the following form in cylindrical coordinates: 

2 = 0; z = l; p cos 9 = 0; psinq) = 0. 

Since the region D is not closed along p, together with the 
boundary conditions (3) it is necessary to take into account also 
the behavior of the field at infinity. Therefore, the potential of 
the field U must be selected in the form of the structure (7) in 
the present case, thus ensuring the requisite behavior of \the poten-.j 
tial and its partial derivatives when p-^<» . The quantities occurring 
in (7-0 have the following form:. 
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V'iei>j+ Vjoii 


(D = (i)i(i)2l 


<i >2 + W| +(l> *■ 

<ai [fl/2Ao/3l‘. <*>2 = — [/ 4 / 5 A 0 M ; 

( (Bo^ [/rAo/aAo/gJ; /i = p — 


fi=ri+hi — p^0: /3=d® — z2^0; 
fi = p — Ta^O; /s = '’2 + ^2 — P^O*. 

fs — (2 — Z + rfa) (Z+rfa — z) ^0; 


/;= j (Z — 2)2^0; /s = p cos (f^O; 
fg = P sin cf^O. 


( 22 ) 


We shall represent <l> In the form of the following approximation: 

p 

O \ ^ % r\ ^ k 

(23) 


0= y; c*„p*(-g-]p, (<p)p,('2), 

where P*, Pi,/*, is some complete system of functions 


Substituting (23) into (7) and performing the operations de- 
scribed above, we dbtaln for finding Cm, a system of the type: 


/lO 


> n 2 I 
-• • .• 


^*0 3 3 V~3T “5T- - 57 -j 

*^/.0 0 0 0 

= ?T i ^ ^ 


(24) 


000 


where Vo and Wma have the form (10a), and in the present case; 


x«»=p*i§-)p/ (<p) Po (^)- 


(25) 


The numerical realization of the system (24), as well as (15) 
also, was accomplished with the help of an electronic computer. 


As an example of the numerical realization of the method de- 
scribed, the structure was selected which is represented in Figure 
la. The solution of the corresponding algorithm was accomplished on 
the electronic computer "Mlnsk-32" for various values of the dimen- 
sions of the bars and the distances between them. The system (15) 
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was solved and the values of the coefficients ■ C*/ (in the two- 
dimensional case) obtained. Then the patterns of the distribution 
of the equlpotentlal surfaces in the region D were constructed in 
accordance with formula (^), taking , into account (8). 




The polynomials of the form: 

X*/=.v'V. {*+7 = 0; 

were chosen as a complete system of functions. 

The set of coordinates functions were chosen subject to the 

condition that their subsequent increase would provide a correction 
to the result not exceeding the prescribed precision of the numeri- 
cal realization of: the potential at any pdlhfe .of the domain D under 
consideration. 


/II 
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The results of the calculation in the form of patterns display- 
ing the distribution of the equlpotentlal surfaces are presented in 
Figure 2a, b, c. The values of the potentials at the corresponding 
boundaries are 1^1 = — 100 . V'2= + I 00 , 1/3= — lOOB, 1/4=0 . 

The equipotentials in Figure 2a, b, c are denoted by numerals, 
representing 7v^~tag e~s" as follows : 

i — 100 ; 2 — 80 ; 2 — 60 ; 4 — 40 ; 5 — 20 ; 6 — 0 ; 7 — 20 ; 5 — 40 ; 9 - 60 ; 

79 — 80 ; 77 — 100 V 

The relative sizes of the regions in Figure 2a: 

«, = 0 , 225 , 02=0,2/5, 03=0.725.04=0.775. 05= 1 ; 6, = 0 , 5 , 62=0.501, 63=0,5. 
64=0.501.65=0.5.66=0.501.67=1. 

The relative sizes of the regions in Figure 2b: 

c, = 0 , 125 , 02 = 0 . 375 , 03 = 0 . 625 , 04 = 0 . 875 , Os=l; 6, = 0 , 1 , 62=0.35, 63=0,1. 

64=0.35, 65=0.1. 66 = 0 . 35 , 67=1. 

The relative sizes of the regions in Figure 2c: 

Oi= 0 , 25 , 02 = 0 , 4 , 03 = 0 , 6 , 04 = 0 , 75 , 05=!; 6i = 0 . 55 , 60 = 0 . 75 , 63 = 0 . 25 , 64 = 

= 0.750,65=0.25,66 = 0.45,67 = 1. 

In conclusion it should be noted that the proposed method, and 
the machine algorithm based on it make it possible to calculate the 
field in structures with elements having arbitrary profiles, and it 
can be used effectively during the development of electron optical 
systems, the focusing devices of superhlgh [frequency generators 
(in particular GDI*), mass spectrometers, and linear accelerators. 

* 

Translator’s note. Expansion unknown. 
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